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ABSTRACT 


The greatest of constant false alarm rate processor (GO CFAR) is a useful 
architecture for adaptively setting a radar detection threshold in the presence of clutter 
edges. The GO CFAR input is often the envelope detected in-phase (1) and quadrature 
(Q) channels of the baseband signal (x, = sqrt (I° + Q*)). This envelope detection can also 
be approximated using x=alI|+b'Q| which requires less complex hardware (a and b are 
simple multiplying coefficients). The envelope GO CFAR processor and several 
envelope approximation GO CFAR processors are compared in terms of the probability 
of detection (PD) performance. Closed-form expressions which describe the PD 


performance are given and their accuracy evaluated. 
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I. INTRODUCTION 


A. GO CFAR PROCESSOR 


The greatest of constant false alarm rate (GO CFAR) processor is a useful 
architecture for adaptively setting a radar detection threshold in the presence of clutter 
edges. The inputs to the processor often involve the envelope detected in-phase (I) and 


quadrature (Q) channels of the baseband signal x, = ii +Q° . This envelope detection 
can be approximated using x=amax{|/|,|Q|}+5min{|/{, |O|} which requires less 
complex hardware and where a and b are simple multiplying coefficients (Pace, 1994 
and Hache, 1994). This approximation can be simplified further to x = alJ| + b|O| which 
reduces hardware requirements even more That is, the two channels do not have to be 


compared to determine the minimum and the maximum. Ironically, this reduction in 
complexity has increased the difficulty of obtaining closed form analytical solutions to 
the GO CFAR performance in terms of probability of false alarm (PFA) and probability 
of detection (PD). Recent studies have used numerical and Monte Carlo methods to 
detail the performance of this processor in terms of the PFA as a function of the 
threshold multiplier. Curve fit solutions have also provided closed form expressions for 
the PFA (Pace, 1994). Closed form expressions for the detection performance however, 


have not yet been examined. 


B. PRINCIPLE CONTRIBUTIONS 


This paper extends the results from previous studies and develops closed form 
expressions for the PD performance as a function of the signal-to-noise ratio (SNR). 
These results provide system designers with an accurate estimate of system performance 
for several detector approximations without requiring numerical analysis or Monte Carlo 
simulations. Since no closed form expressions for the PD exist, Monte Carlo simulations 
are first used to obtain a full set of the PD curves using a PFA of 10“ for the number of 


reference cells n=1,2.4,8.16.32 and for various values of multiplying coefficients a and b. 











Various curve-fit techniques are then used to obtain closed form expressions for 
the PD vs. SNR curves. Several formulations are investigated. The most promising 
expression involves the use of the erf function which produces a detection performance 
within +0.015 of the actual PD. 

Next, the coefficients for each envelope approximation PD performance are 
plotted as a function of the number of cells n. This allows the extraction of the 
appropriate coefficients for any number of reference cells n between 0 and 32. This in 


turn provides a quick and relatively accurate expression for the PD as a function of the 


SNR for the various detection approximations a and b.. 


C. THESIS OUTLINE 


Chapter II starts with a review of basic signal detection techniques, an 
introduction to the GO CFAR processor, and the concepts of noise and probability. It 
continues with a description of envelope detection as well as several approximations to 
the envelope detector and concludes with an analysis of the GO CFAR processor 
performance in terms of the PFA and the PD. Chapter Ill provides a discussion of the 
Monte Carlo techniques and the simulation results. Chapters I'V presents the derivation of 
closed form expressions for the detection performance and provides tables of the 
corresponding coefficients. Chapter V presents the plots of the closed form expression 
coefficients as a function of the number of cells n. The appendix contains the code 


developed for this research. 











II. CFAR PROCESSOR 


A. BASIC SIGNAL DETECTION PRINCIPLES 


The detection process for a typical radar system involves an inherently 
straightforward comparison between the received signal and a reference value. In a 
simple detector, the reference value is a fixed threshold which remains unchanged 
through the detection process. If the received signal is above the threshold, a target is 
assumed to be present, otherwise the signal is ignored. Not surprisingly, this is referred to 


as threshold detection. Figure 1 illustrates this concept. 


Threshold level 





Figure |. 
Typical Envelope of a Radar Receiver Output as a Function of Time. (After Skolnik, 1980) 


A low threshold increases the chance that noise alone will rise above the threshold and 
generate false alarms. On the other hand, a high threshold will result in missed targets. 
Consequently, the selection of a proper threshold is a compromise between reducing the 
PFA and increasing the PD. 

To make things more difficult, subtle changes in the threshold can drastically 
change the overall system PFA. Although these changes can be manually compensated 
by trained operators (indeed this is how old systems operated), current and future 
processors require near-instantaneous changes in threshold levels to maintain a constant 
PFA to prevent overloading automatic detection and tracking systems (Skolnik, 1980). 


One such processor is the GO CFAR or Greatest-Of, Constant False Alarm Rate device. 
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B. GOCFAR PROCESSOR 


In the GO CFAR processor, the background noise or clutter in the vicinity of a 
target cell is sampled, then the threshold is adjusted to provide a constant false alarm 


rate. Figure 2 depicts a how a typical radar might view a target in a particular region in 


space. 


Range cells ————__—______—-_> 





Figure 2. Typical Radar View of Target Region. 


For sake of simplicity, assume that a target is located completely within a single range 
cell, while the remaining cells contain only thermal noise. Once received, the signal is 
frequently manipulated to create a quadrature channel which is 90° out of phase with the 
in-phase signal. The purpose of the quadrature channel is to retain the phase information 
and eliminate the effects of blind phases (Skolnik, 1980). 

A schematic diagram of the GO CFAR processor is shown in Figure 3. The input 
1 / Q thermal noise samples are considered to be normally distributed N(0,1) and 


detected as either x= //?+Q? or x=all| +b|Q| in the test and reference cells. Both 
reference cell neighborhoods contain n cells which are used to determine the noise power 


levels y,, y,. The detector threshold voltage V, is obtained by choosing the greater of y, 
and y,, normalizing by the number of reference cells n, and multiplying by the threshold 


multiplier 7. A target is declared if the amplitude of the test cell is greater than /’,. 
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Figure 3. GO CFAR Processor (After Skolnik, 1980) 





The threshold multiplier (T) is normally set a priori to achieve a desired false 
alarm rate (Pace, 1994). For any fixed T, increasing the number of cells sampled (7) 
corresponds with a decrease in the false alarm rate. The probability of detection, on the 
other hand, is based upon the SNR of the received signal and is normally outside the 
realm of the designer's control. | 

For mathematical convenience, CFAR processors are frequently analyzed using 
the envelope detector x- = {[P+@ to detect the In-phase (1) and Quadrature (Q) 


channels of the baseband signal. To reduce hardware complexities created by the use of 


squares and square roots in the envelope detector, numerous envelope approximations 
have been developed in the form x =al/| +b|O| where a and b are simple multiplying 
coefficients (Pace, 1994). These simplifications result in a slight loss in sensitivity, a 


difference in PFA performance (Pace, 1994) and consequently, a decrease in detection 
probability. This study analyzes the detector performance of the GO CFAR processor 
using the envelope detector and several envelope detector approximations. Since 


numerical integration of the PD expressions require considerable time, Monte Carlo 








simulations are used to obtain the results. Closed form expressions for the PD are then 
found using curve fitting techniques. 
C. NOISE AND PROBABILITY 

Radar signal analysis requires an initial understanding of noise which is a random 
phenomenon. 

1. Probability Density Functions 


Many sources of noise can be modeled by a gaussian or normal density function, 


also commonly referred to as the bell curve as illustrated below in Figure 4, 





Figure 4. Gaussian, Normal or Bell Curve 





and given by 
| _Gen)? 
x)=—=—e 2° 3.1 
AX) = ie (3.1) 


where 1=mean and o = standard deviation, N(,0). The x-axis represents signal 
amplitude, the y-axis the probability of a signal associated with that amplitude Using 


this model, noise is comprised primarily of low amplitude noise centered around zero, 











whereas the probability of higher amplitude noise drops off exponentially with 
increasing amplitude. 
2. Cumulative Distribution Functions 


The cumulative distribution function (cdf) is the probability that the value x 1s 


less than some specified value and is given mathematically as the integral of the density 


function 
Fy) = ie Iix)dx. (3.2) 


By convention, cdfs are represented by capital letters, while pdfs are represented by 


lower-case letters. The cdf of the density function in Figure 4 appears in Figure 5. 





Figure 5. Typical Cumulative Distribution Function. 








3. Summation of Density Functions 


Given two different functions of similar, but independent random variables the 
pdf of the output is equal to the convolution of the density functions (Papoulis, 1984): 


f-(z) = | __ flz—yfy(v) dy (3.3) 


Note that if f,(x) = 0 for x <0 and f,(0) = 0 for y <0, then f,(0) = 0 forz <0 and 


f(z) = | f(z —yfyo)dy z>0. (3.4) 


4. Maximum of Two Density Function 


Given two functions of independent random variables f,(x) and f,(v), cdf of the 


maximum of the two functions is given as (Papoulis, 1984): 
F(z) = F(z) Fy) (3.5) 
and the corresponding pdf is given as: 


f(z) =f (2) F(z) + fp) Fx (2). (3.6) 
If the pdfs of f,(x) and f,(v) are equivalent, then the pdf simplifies to 


f(z) = 2fx(2)F x(z) = 2fv@)F'y(Z). (3.7) 
5, Linear transformation of density functions 


For the linear transformation y=ax where x is a random variable with pdf f,(x), 
the pdf of y, f(y) is (Papoulis, 1984): 


fv) = rafx(a). 38) 











6. Absolute value of a density function 


For the absolute value of a function of a random variable y=|x|, the resulting pdf 


is given as 


hO) =A) +h). (3.9) 


The relationships discussed above are needed to determine the density function for the 


detector approximation being examined. 


D. ENVELOPE DETECTOR AND APPROXIMATIONS 


The input to the GO CFAR processor is composed of a combination of in-phase 


and quadrature signals which are 90° out of phase and of the form /=Acos()+ x and 
O=Asin(o)+y where A is the signal amplitude, o is the signal phase and x, v are 


Gaussian random variables with zero mean and variance of one N(0,1). The associated 
pdfs of the absolute value of these two signals and associated curves, with o=7/8 are 


shown below: 


_(x-Acos (gy)? _ (-¥-A.c0s(6))? 
2 +e 2 
fin(x, A, 0) = ————>__ «> 0 (3.10) 


(20 





Figure 6. PDF of |I| with o=1/8 








_(eAsin(g)? —_ (cxe-dsin(g))? 
2 


fig\(x,4,o) = —————__ x>0 G.11) 


f(x, @) 0.4 





Figure 7. PDF of |Q| with o=p/8 


Using Equation (3.8) to factor in the a and b coefficients for the envelope approximation, 
we obtain the following pdfs: 


( ~_Acos (6) | ; 


fall (% 0) == io x >Q (3.12) 


la| 


(2-4 cos () : 





_e. 2 +e. 2 
fojo|; 9) = mS x>0 (3.13) 





As previously discussed, the envelope detector is given by x- = (P. + Q* while 
the approximation formula is x = a|/| + b|Q|. Seven sets of multiplying coefficients for a 


and b which provide reasonable approximations to the envelope detector are given in 


Table 1 below. 
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Table 1. 
Multiplying Coefficients a and b for Seven Envelope Approximations (After Filip, 1976) 


E. PROBABILITY OF FALSE ALARM VS THRESHOLD MULTIPLIER 


Using the results of previous sections, the PFA of the GO CFAR processor 1s 
shown. The GO CFAR processor is reproduced in the figure below. 





Figure 8. GO CFAR (After Skolnik, 1980) 


Since we're interested in the PFA, we start with the assumption that no target 1s 
present in the test cell and ask ourselves: What is the probability that the signal Tz/n is 
higher than the test cell output x? The answer to this question gives us the probability of 
false alarm. We start with the pdf of the input signal, step through the GO CFAR model 
calculating corresponding pdfs at each junction, and end up with the PFA at the output of 


the GO CFAR. 


1] 





With no target present (A=0), the pdfs of |I| and |Q| are greatly simplified. In fact, 


they are the same. 


fin(x) =f\o| = Be | | (3.14) 


For the envelope approximation x = al/| +b|Q|, we use our previous results to 
calculate the pdf of the input as the convolution of the pdfs of a|I| and b|Qj. The resultant 





pdf becomes 
ee 
2 2 a2+6? | a x b x 
fx(x) = e | erf{| ——_=—_- |+ er (3.15) 
x(X) 1(a*+b*)  [2a2+b?) * [2(a?+b?) 
where erf{x) is the error function given by 
Lop 
erf(x)=—J|,e 2 dy. 3.16 
Ax) = = Je 7 dy (3.16) 


Assuming the input contains only noise with no target present, Equation (3.16) also 
represents the pdf of a test cell output. 

The summation of input signals which result in the outputs yl and y2 are a 
function of the number of reference cells n. Using the results of Equation (3.4). we note 
that f.(yl) and f,,(y2) are both n-fold convolutions of f(x) for which there is no 


closed-form solution. 
The pdf of z after the Greater-Of block can be determined using Equation (3.7): 


fz(2) = 2fvr (2) 1 (2) = 2fr2(2)F 22) (3.17) 


where F(z) is the cumulative distribution function of f,(z). The PFA at the comparator 


output is then 
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PFA(T) = | flat | "fda faz. 3.18) 
Substituting in Equation 18 and integrating F,, we get 


PFA(T) =2 f° fa {J fuel] [ K@dr |de 3.19 


where f(A) is given Equation (3.15). Although this equation has no closed form solution, 
it can be calculated numerically. Also, closed form expressions for the PFA are given in 
Pace. These results are used to extract the corresponding threshold multipliers used in the 


Monte Carlo simulations for the detection performance. 


F. PROBABILITY OF DETECTION VS. SNR 


The determination of the probability of detection follows a similar line of 
analysis as in the previous section. We now assume however, a target is present in the 
test cell. In addition, the threshold multiplier T is no longer a variable, but a constant 
which corresponds with a particular PFA. This time, we ask the question: What is the 
probability that the signal from the test cell containing the target is higher than the signal 
T2/n appearing at the comparator? The answer to this question gives us the probability of 
detection. 

Our first step is to calculate the pdf of the input signal x using the envelope 
approximation. Using Equation (3.4) once again, we can calculate the pdf of the 
envelope approximation as the convolution of the pdfs of all| and 5|Q| with a target 


present (Equations (3.12) and (3.13) respectively): 


(+ 4008(6))” (24100816) (a4 sino) . (4-1 sin(o)) 
Ae. 2 +e 2 e_ : +e_ : 
—{ @—* te =n SS 
fx(x, A, O) : al (an \o| (2m - 


(3.20) 
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Substituting in the following relation of SNR to signal amplitude A (Wilson, 1982), 


SNR = 4 


and performing the integration. we obtain the following closed-form solution. 


_(B-x)? 

OF Io Et) ef Ss) 
cs pet) haa) 
() E+H-I 
(-S") E+H-I 

ef a) {et) 


with the following substitutions: 











f(x, SNR, >) = A 








= /2SNR (acos() + b sin(9)) 
C = J/2SNR (acos() — bsin()) 
D=a’*+b? 

= /2SNR b*cos($) 
F = /2SNR a’sin (4) 
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(3.21) 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


(3.28) 











G = J2SNR ab cos(6) (3.29) 
H = J2SNR absin(0) (3.30) 
[=ax (3.31) 
J = bx (3.32) 


The PD equation is similar in form to the PFA equation, Equation (3.19), but T is 


now constant and derived for a specific PFA. 


PD(SNR,4)= J 2fnl{f- fr ©rae}] [" O.SNR, O)dA [az 2.33) 


Since the phase of the target signal is also random, the PD must be averaged over the 


phase as 
PD(SNR) =f 2 fala fudeBl] J” AO, SNR,AVMA |dedo. 3.34 


Once again, there is no closed form expression for the PD equation, but it can be 
evaluated using numerical integration or Monte Carlo simulations. The latter approach 1s 
the subject of the following chapter. This is necessary in order to develop the closed form 


expressions to approximate this detection performance. 











Ill. MONTE CARLO SIMULATIONS USING MATLAB 


A. MATLAB PROGRAM 


Monte Carlo simulations are conducted using the MATLAB code located in the 
appendix. The goal is to obtain the PD vs. SNR curves associated with a PFA value of 
10“ for each envelope approximation. The PFA chosen is a value generally reported in 
the literature for CFAR comparisons. For both the envelope approximation and envelope 
detector, the following number of reference cells n =1,2,4,8,16,32 are used. Threshold 
multipliers (7) for each particular envelope approximation, the number of reference cells 
(n) and the desired PFA, are obtained by interpolating the PFA vs. the threshold 
multiplier data from a previous study (Pace, 1994). The tables below summarize the 


threshold multipliers used for each Monte Carlo simulation. 





Table 2. Threshold Values for Envelope Approximation a=1, b=1 


Table 3. Threshold Values for Envelope Approximation a=1, b=1/2 
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Table 7. Threshold Values for Envelope Approximation a=0.948, b=0.393 
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Table 9. Threshold Values for Envelope Detector {P +Q 


B. PROGRAMMING TECHNIQUES 


To simulate noise-only inputs to the reference cells, the program creates matrices 
of random numbers N(0,1) with dimensions (n x N), where n = number of reference cells 
per side and N = number of simulations desired. The absolute values of these matrices 
are summed across the n variable to create a vector of dimension (1 x N), then multiplied 
by the appropriate coefficients a or b and added together. This creates two vectors of 
dimension (1 x N) which represent the summation of reference cells on the left and nght 
side of the test cell. These two vectors are compared, and for each N, the greater of the 
two cells is placed into a new vector z of dimension (1 x N). To create the comparison 
voltage Vt we divide by the number of reference cells n, and multiply by the threshold 
multiplier T. The test cell signal, Vtest, is created in a similar manner, but this time we 
add a signal with amplitude A which corresponds to the SNR to be tested. 

The next step compares the two vectors Vt and Vtest and calculates the 


percentage of times that Vtest > Vt. For each SNR tested, the process is repeated for 
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0<<n/2 since the phase of noise is equally distributed. The result is a probability of 
detection for a single SNR. This must be repeated for the entire range of SNRs to be 


tested. 

After a number of trial runs, it was determined that 400,000 simulations per SNR 
with an interval of 0.5 dB provided reasonably smooth curves for our analysis. We also 
took advantage of the fact that at some point, when the SNR becomes large enough, the 
PD reaches 1.0. Thus, the algorithm breaks the loop after five successive PDs of 1.0 to 


remove unnecessary computational loops. The PD plots are given below. 
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Figure 9: PD vs. SNR for a=1, b=1, PFA=10” 
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Figure 11. PD vs. SNR using Envelope Approximation, PFA=10™, a=1, b=1/4 
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Figure 13. PD vs. SNR using Envelope Approximation, PFA=10™, a=31/32, b=3/8 
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Figure 15. PD vs. SNR using Envelope Approximation, PFA=10", a=0.96043, b=0.39782 
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Figure 16. PD vs. SNR using Envelope Detector, PFA=10" 
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IV. APPROXIMATION FORMULAS FOR PD 


This chapter outlines the derivation of closed-form, curve fit expressions for the 


various PD vs. SNR curves. 


A. METHODS OF CURVE FITTING 


Numerous methods were attempted to curve fit the Monte Carlo results and 
included the use of tanh, third-order exponentials, and erf functions. The most promising 


involved the use of the ef function. 


B. ERF FUNCTION APPROXIMATION 


The erf or error function, previously given in Equation (3.16), is repeated here 


along with a plot of the function for convenience 


2 
ee 
erfix) = == J) ¢ 2 dy (4.1) 





Figure 17: Plot of Erf Function 
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The erf function, due to its similar appearance to the PD curves in the previous section, 1s 


a natural candidate as a curve fit solution. By rewriting the function in the form 


Kx) = 0.5( enf =] za | (4.2) 


we can control both the curvature and axes with the proper values of c/ and c2. Using the 
PD vs. SNR curve for n=1 as an example, we can obtain a reasonably good fit using 


=22.5 and d=5.26. 
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Figure 18: Comparison Between Actual PD Curve and Curve Fit Solution Using Evf Function 


Despite the symmetric nature of the erf function, we are still able to produce a reasonable 
curve fit to the asymmetric solid PD curve (the line we are attempting to curve fit) and, 


at the same time, obtain a relatively low residual across the entire SNR axis. The residual 


is defined as 
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Residual(SNR) = |PD(SNR) — f(SNR)| (4.3) 


and describes the error or difference between the two curves. The residual plot is shown 
below. 
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Figure 19: Residual of Erf Function Approximation 
In this case, the maximum residual equals 0.0218. That is, for any given SNR, the erf 
approximation yields a PD result within 0.0218 of the actual value. 


To improve this result, the curve-fit solution is split into two parts: PDs > 0.5 and 
PDs < 0.5. A curve fit solution is obtained for the upper portion of the curve above a PD 


of 0.5, and a second solution obtained for the /ower curve. The formulas are given as: 


PD<0.5 PD = 0.5(erf(=5) + 1) (4.4) 


PD > 0.5 PD = 0.5(erf( =) +1) (4.5) 
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where the coefficients c/, c2, c3 are determined experimentally. 


In this case, the maximum error is reduced to 0 0116, clearly a better solution. 
Plots of the improved solution and corresponding residual curve for a=1, b=1, n=1, are 


given in the following figures. 
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Figure 21: Residual of Approximation 
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c. TABLES OF COEFFICIENTS 


The following tables provide the associated coefficients for each of the seven 


envelope approximations and envelope detector. 





Table 12. Coefficients for PD Curve PFA=10", a=1, b=1/4 
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Table 13. Coefficients for PD Curve PFA=10", a=1, b=3/8 
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Table 14. Coefficients for PD Curve PFA=10", a=3 1/32, b=3/8 










Table 16: Coefficients for PD Curve PFA=10", a=0.96043, b=0.39782 
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Table 17. Coefficients for PD Curve PFA=10%, Envelope Detector 
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V. PD APPROXIMATION COEFFICIENTS 


A curve fit of the coefficients as a function of the number of reference cells n, 
allows system designers to extract the appropriate coefficients and obtain approximate 
PD vs. SNR curves for any number of reference cells n. The complete set of coefficient 


plots are given below. 
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Figure 22: Coefficient C1 as a Function of n 
For Envelope Approximation a=1, b=1 
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Figure 23: Coefficients C2 (Solid Line) and C3 ( Dashed Line) as a Function of n 
For Envelope Approximation a=1, b=1 
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Figure 24: Coefficient C1 as a Function of n 
For Envelope Approximation a=1, b=1/2 
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Figure 25. Coefficients C2 (Solid Line) and C3 (Dashed Line) as a Function of 
For Envelope Approximation: a=1, b=1/2 
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Figure 26: Coefficient C1 as a Function of n 
For Envelope Approximation a=1, b=1/4 
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Figure 27. Coefficients C2 (Solid Line) and C3 (Dashed Line) as a Function of n 
For Envelope Approximation: a=1, b=1/4 
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Figure 28: Coefficient C1 as a Function of n 
For Envelope Approximation a=1, b=3/8 





3 

2 

2 

og 

20 25 
Number of test cells (n) 
Figure 29. Coefficients C2 (Sc >: and C3 (Dashed Line) as a Function of n 
For Envelop. . . -~7ximation: a=1, b=3/8 
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Figure 30: Coefficient Cl as a Function of n 
For Envelope Approximation a=3 1/32, b=3/8 
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Figure 31. Coefficients C2 (Solid Line) and C3 (Dashed Line) as a Function of n 
For Envelope Approximation: a=31/32, b=3/8 
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Figure 32: Coefficient C1 as a Function of n 
For Envelope Approximation a=0.948, b=393 
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Figure 33: Coefficients C2 (Solid Line) and C3 (Dashed Line) as a Function of 
For Envelope Approximation: a=0.948, b=0.393 
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Figure 34: Coefficient C1 as a Function of n 
For Envelope Approximation a=0.96043, b=0.39782 
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Figure 35. Coefficients C2 (Solid Line) and C3 (Dashed Line) as a Function of n 
For Envelope Approximation: a=0.96043, b=0.39782 
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Figure 36: Coefficient C1 as a Function of n 
For Envelope Detector 


Coeflicients (c2 and ¢3) 


Number of test cells (n) 





Figure 37. Coefficients C2 (Solid Line) and C3 ( Dashed Line) as a Function of n 
For Envelope Detector 
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VI. CONCLUSIONS 


Monte Carlo simulations yield accurate PD vs. SNR curves for envelope 
detection approximation GO CFAR processors. These curves can be closely 


approximated with closed form expressions using the erf function in the form: 


po<0s =PD=0.5(erf(= 5") +1) (6.1) 
po>0s PD=0.5(erf(@S*)+1) (6.2) 


where the coefficients cl, c2 and c3 are chosen for a specific envelope approximation 
and number of reference cells n. These PD approximations compare favorably with the 
results of Monte Carlo simulations; the maximum residual in all cases is less than 0.028. 
A curve fit of the coefficients cl, c2 and c3 as a function of the number of 
reference cells provides a means of quickly extracting the coefficients for any number of 
reference cells 1<n<32 for a PFA=10~. In conjunction with the equations above, these 


coefficients produce accurate PD curves as a function of SNR. 
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APPENDIX - MATLAB PROGRAMS 


Programs to Run Monte Carlo Simulation 
runmont.m _ Batch file to run multiple simulations 


setvar.m Program to set variable for each loop 
montesim.m Main program to run simulations of envelope approximations 


monteenv.m Main program to run simulation of envelope detector 


Programs to run Plots of Data 
plotmont.m Plots data from *.dat files 


getdata.m Gets data to plot graphs 
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RUNMONT.M 


% FILE NAME: RUNMONT.M 


% Batch file to run all Monte Carlo Simulations for PFA=1e-4 


clear 

% Input parameters 

% a: coefficient a 

% b: coefficient b 

%o n: number of cells per side 


% PFA: Required PFA 

% T: Threshold 

% N: number of monte carlo simulations per inner loop 
% MAXSNR: Maximum SNR in dB 


Ose oten a eee ee ae 
% Set variables using data from pfa_a*.dat at PFA=le-4 
% SNR/PD data saved to file MC1_4.dat 

set_var; % Clear all variables and reset to clear memory 
a=1; b=1; 

nTABLE={ 1.000 2.000 4.000 8.000 16.00 32.000); 
TTABLE=[12.287 5.952 4.423 3.891 3.680 3.5975]; 
montesim; 

save MC1_4.dat mctemp -ascil 


% Set variables using data from pfa_b*.dat at PFA=le-4 
% SNR/PD data saved to file MC2_4.dat 

set_var; % Clear all variables and reset to clear memory 
a=1; b=0.5; 

nTABLE=[ 1.000 2.000 4.000 8.0000 16.00 32.000]; 
TTABLE=[13.430 6.515 4.773 4.1525 3.903 3.8040]; 
montesim; 

save MC2_4.dat mctemp -ascii 

Ee EE SOM Re RE Ce EEN ete eee 
% Set variables using data from pfa_c*.dat at PFA=le-4 
% SNR/PD data saved to file MC3_4.dat 

set_var; % Clear all variables and reset to clear memory 
a=1; b=0.25; 

nTABLE=[ 1.000 2.000 4.000 8.0000 16.00 32.000]; 
TTABLE=[16.685 7.915 5.530 4.6750 4.335 4.197]; 
montesim; 

save MC3_4.dat mctemp -ascii 

| LE ae SP APE Ree cea eee ee ee er eee 
% Set variables using data from pfa_d*.dat at PFA=le-4 
% SNR/PD data saved to file MC4_4.dat 

set_var; % Clear all variables and reset to clear memory 
a=1; b=3/8; 

nTABLE={ 1.000 2.000 4.000 8.0000 16.00 32.000]; 
TTABLE=[14.530 7.018 5.064 4.3610 4.076 3.9636]; 
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montesim; 
save MC4_4.dat mctemp -ascii 


% Set variables using data from pfa_e*.dat at PFA=le-4 
% SNR/PD dita saved to file MCS_4.dat 

set_var; % Clear all variables and reset to clear memory 
a=3 1/32; b=3/8; 

nTABLE=[ 1.000 2.000 4.000 8.0000 16.00 32.000]; 
TTABLE=[14.399 6.960 5.032 4.3380 4.058 3.9465]; 
montesim; 

save MC5_4.dat mctemp -ascil 


% Set variables using data from pfa_f*.dat at PFA=le-4 
% SNR/PD data saved to file MC6_4.dat 

set_var; % Clear all variables and reset to clear memory 
a=0.948; b=0.393; 

nTABLE={ 1.000 2.000 4.000 8.0000 16.00 32.00]; 
TTABLE={14.112 6.835 4.960 4.2870 4.017 3.907]; 
montesim; 

save MC6_4.dat mctemp -ascil 


% Set variables using data from pfa_g*.dat at PFA=le-4 
% SNR/PD data saved to file MC7_4.dat 

set_var; % Clear all variables and reset to clear memory 
a=0.96043; b=0.39782; 

nTABLE=[ 1.000 2.000 4.000 8.0000 16.00 32.00]; 
TTABLE=[14.114 6.834 4.960 4.2880 4.017 3.908}; 
montesim; | 
save MC7_4.dat mctemp -ascii 

Of anmnnennnrnnmnersceceerannnnnnnennnennesennnnceenanawanns 
% Set variables using data from pfa_r*.dat at PFA=le-4 
% SNR/PD data saved to file MCr_4.dat 

set_var; % Clear all variables and reset to clear memory 
nTABLE=[ 1.000 2.000 4.000 8.0000 16.000 32.000); 
TTABLE={ 11.962 5.751 4.251 3.7290 3.5215 3.4395]; 
monteenv; 

save MCr_4.dat mctemp -ascil 
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SET_VAR.M 


% FILE NAME: SET_VAR.M 
% File to set variables 


clear global 

clear all 

% Input parameters 

% a: coefficient a 

% b: coefficient b 

% n: number of cells per side 

% PFA: Required PFA 

% T: Threshold 

% Nmax: Max number of monte carlo simulations per angle phi 
% MAXSNR: Maximum SNR in dB 


angle=pi/2;% Set max angle to sample (0 to pi/2) 

phi=0:angle/40:angle;% Set phi to vector of angles from 0 to pi/2, 40 steps 
PFA=le-4; 

Nmax=10000; 

MAXSNR=40; 
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MONTESIM.M 


% FILE NAME: MONTESIM.M 
% Generic Monte Carlo Simulation 


% 

“% Input parameters 

% a: coefficient a 

% b: coefficient b 

% n: number of cells per side 


% PFA: Required PFA 
% T: Threshold 


% Nmax: Max number of monte carlo simulations per angle phi 

% N: Size of random matrix (nxN)=Nmax | 

% MAXSNR: Maximum SNR in dB 

% nTABLE: Values of n=number of cells per side 

% TTABLE: Values of T=corresponding threshold setting for given n and PFA 
% angle: Max angle to sample (0 to pi/2) 

% phi: vector of angles from 0 to pi/2, 40 steps 


for num=1:length(nTABLE), 

skip=0; % Sets flag to zero (Used to remove unnecessary loops) 

count=0; % Sets count to zero (Used to remove unnecessary loops) 
n=nTABLE(num); 

T=TTABLE(num); 

N=round(Nmax/n); % Size of random matrix (nxN)=Nmax 


SNRdB=0:0.5:MAXSNR;% Set Signal-to-Noise ratio in dBs (increments of 0.5 dB) 
PD=zeros(1,length(SNRdB)); 


i ia Sa as eae aan a eenad 

for i=1:length(SNRdB); % Performs loop over all SNRs 

if skip==1 
break 

end 

for j=1:length(phi), % Performs loop over all angles phi: 0 to pi/4 (20 steps) 
SNR=10*(SNRdB(i)/10);% Converts SNR from dBs to linear scale [1,1] 
A=sqrt(2*SNR); % Converts SNR value to required A [1,1] 
In]=randn(n,N); % 1 on left side, noise only [nxN] 
In2=randn(n,N); % I on right side,noise only [nxN] 
Qn1=randn(n,N); % Q on left side, noise only [nxN] 
Qn2=randn(n,N); % Q on left side, noise only [nxN] 

if n=] 


y 1=a*abs(In])+b*abs(Qn1); 

y2=a*abs(In2)+b*abs(Qn2); 

else 

y1=a*sum(abs(In1))+b*sum(abs(Qn1));  % Sum of left side noise inputs [1xN] 
y2=a*sum(abs(In2))+b*sum(abs(Qn2)); % Sum of right side noise inputs [1xN] 
end 
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z=max(y1,y2);% GO output [1xN] 

Vt=T *z/n; % Comparison voltage [1xN] 
A ee aT eS 

% Test voltage in test cell with signal [1,N] 

% z=a*|I|+b*|Q| or z=a*|Acos(phi)+x|+b*|Asin(phi)t+y| 

Vtest=a*abs( A *cos(phi(j))+randn(1,N))+b*abs(A *sin(phi(j))+randn(1,N)); 
ER ee nee re orn 

% Calculate percentage 

% Add up all percentages for each value of phi 
PD(i)=PD(i)+length(find((Vtest-Vt)>0))/(N*length(phi)); % [1 x # of SNRs} 
end | 

% deletes unnecessary loops if PD=1 for greater SNRs 

if PD(i)>=1 

count=count+1; 

if count>5 

skip=1; 

end 

end 
end 





Wqn--= Break branches to here 


for k=1:length(SNRdB), 
if k>=i 

PD(k)=1; 

end 


end 
% Store SNR/PD data to file named mctemp (temporarily) 


mctemp((2*num)-1:2*num,:)=[SNRdB;PD]; 
end 
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MONTEENV.M 


% FILE NAME: MONTEENV.M 
% Generic Monte Carlo Simulation for envelope detector 


% 
“% Input parameters 
% n: number of cells per side 


% PFA: Required PFA 
% T: Threshold 


% Nmax: Max number of monte carlo simulations per angle phi 
% N: Size of random noise matrix (nxN)=Nmax 

% MAXSNR: Maximum SNR in dB 

% nTABLE: Values of n=number of cells per side 


% TTABLE: Values of T=corresponding threshold setting for given n and PFA 
% angle: Max angle to sample (0 to pi/2) 
% phi: vector of angles from 0 to pi/2, 40 steps 


for num=!1 :length(nTABLE); 

skip=0; % Sets flag to zero 

count=0;% Sets count to zero 

n=nTABLE(num); 

T=TTABLE(num); 

N=round(Nmax/n); % Set size of random noise matrix (nxN)=Nmax 


SNRdB=0:0.5:MAXSNR;% Set Signal-to-Noise ratio in dBs (increments of 0.5 dB) 
PD=zeros(1,length(SNRdB)); 


for i=1:length(SNRdB), % Performs loop over all SNRs 
if skip==1 

break 

end 


for j=1:length(phi), % Performs loop over all angles phi: 0 to pi/4 (20 steps) 


SNR=10“(SNRdB(1)/10); % Converts SNR from dBs to linear scale [1,1] 
A=sqrt(2*SNR); % Converts SNR value to required A [1,1] 
In}=randn(n,N): % I on left side, noise only [nxN] 

In2=randn(n,N); % | on right side,noise only [nxN] 

Qn1=randn(n,N); % Q on left side, noise only [nxN] 

Qn2=randn(n,N); % Q on left side, noise only [nxN] 

if n==] 


y l=sqrt((In1.42)+(Qn1.%2)); 
y2=sqrt((In2.*2)+(Qn2.%2)); 


else 

y l=sum(sqrt((In1.42)+(Qn1.42))); % Sum of left side noise inputs [1xN] 
y2=sum(saqrt((In2.*2)+(Qn2.%2))); % Sum of right side noise inputs [1xN] 
end 


z=max(yl,y2); % GO output [1xN] 
Vt=T *z/n; % Comparison voltage [1xN] 


% Test voltage in test cell with signal [{1,N] 
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% z=sqrt[(1*2)+(Q*%2) or z=sqrt[(Acos(phi)+x)*2+(Asin(phi)+y)*2] 
Vtest=sqrt((A *cos(phi(j))+randn(1,N)).*2+(A*sin(phi(j))+randn(1,N)).”2); 
Of nenennnnnann mann ner n enna nen nen en nnn neemnn 


% Calculate percentage 
% Add up all percentages for each value of phi 
PD(i)=PD(i)Hength(find((Vtest-Vt)>0))/(N*length(phi)); % [1 x # of SNRs] 


end 
% deletes unnecessary loops if PD=1 for greater SNRs 
if PD(i)>=1 
count=count+1; 
if count>5 
skip=1; 
end 
end 
end 


%/,...- Break branches to here 


for k=1:length(SNRdB), 
if k>=1 

PD(k)=1; 

end 


end 
% Store SNR/PD data to file named mctemp (temporarily) 


mctemp((2*num)-1:2*num,:)=[SNRdB;PD}; 
end 
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PLTMONT.M 


% FILE NAME: PLTMONT.M 

% Batch file to plot PD vs. SNR curves from Monte Carlo Simulations 
% Plot PD vs SNR curve from mcl_4.dat 

load mci_4.dat 

save mc.dat mcl_4 -ascil 

var_title="PD vs. SNR a=1,b=1,PFA=le-4, n=1,2,4,8,16,32', 
get_data 

% Plot PD vs SNR curve from mc2_4.dat 

load mc2_4.dat 

save mc.dat mc2_4 -ascil 

var_title='PD vs. SNR a=1,b=1/2,PFA=le-4, n=1,2,4,8,16,32'; 
get_data 

% Plot PD vs SNR curve from mc3_4.dat 

load mc3_4.dat 

save mc.dat mc3_4 -ascil 

var_title='"PD vs. SNR a=1,b=1/4,PFA=le-4, n=1,2,4,8, 16,32’; 
get_data 

% Plot PD vs SNR curve from mc4_4.dat 

load mc4_4.dat 

save mc.dat mc4_4 -ascii 

var_title='"PD vs. SNR a=1,b=3/8,PFA=1e-4, n=1,2,4,8,16,32'; 
get_data 

% Plot PD vs SNR curve from mc5_4.dat 

load mc5_4.dat | 

save mc.dat mc5_4 -ascli 

var_title="PD vs. SNR a=31/32,b=3/8,PFA=1e-4, n=1,2,4,8,16,32'; 
get_data 

% Plot PD vs SNR curve from mc6_4.dat 

load mc6_4.dat 

save mc.dat mc6_4 -ascii 

var_title='PD vs. SNR a=0.948,b=0.393,PFA=1e-4, n=1,2,4,8,16,32'; 
get_data 

% Plot PD vs SNR curve from mc7_4.dat 

load mc7_4.dat 

save mc.dat mc7_4 -ascii 

var_title='PD vs. SNR a=0.96043,b=0.39782,PFA=1e-4, n=1,2,4,8,16,32'; 
get_data 

% Plot PD vs SNR curve from mecr_4.dat 

load mcr_4.dat 

save mc.dat mcr_4 -ascil 

var_title='PD vs. SNR for envelope detector,PFA=le-4, n=1,2,4,8,16,32'; 
get_data 


53 








rear a ee eee ————————————————————— 


54 











GET_DATA.M 


% FILE NAME: GET_DATA.M 
% Batch file to plot PD vs. SNR curves from Monte Carlo Simulations 


% Get data 

clf 

load mc.dat 
num_plot=length(me(:,1))/2; 
hold on 

grid on 

for i=1:num_plot; 
plot(mc((2*i)-1,:),mc((2*1),:)) 
end 

xlabel('Signal-to-Noise Ratio (dB)') 
ylabel('Probability of Detection’) 
title( var_title) 

% print -dwin 

pause 
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